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ABSTRACT 

As part of the VLT/X-shooter science verification, we obtained the first optical medium- 
resolution spectrum of a previously identified bright O-type object in NGC 55, an LMC-like 
galaxy at a distance of ~2.0Mpc. Based on the stellar and nebular spectrum, we investigate 
the nature and evolutionary status of the central object(s) and its influence on the surround- 
ing interstellar medium. We conclude that the source, NGC55C1_31, is a composite object, 
likely a stellar cluster, which contains one or several hot (T e ff - 50000 K) WN stars with 
a high mass-loss rate (~3 x 10~ 5 M yr _1 ) and a helium-rich composition (Nne/NH = 0.8). 
The visual flux is dominated by OB-type (super)giant stars with T^g < 35 000 K, solar he- 
lium abundance (Nue/NH = 0.1), and mass-loss rate ~2 x 1O~ 6 M yr _1 . The surround- 
ing Hn region has an electron density n e < 10 2 cm -3 and an electron temperature r(Oin) 
- 1 1 500 ± 600 K. The oxygen abundance of this region is [O/H] = 8.18 ± 0.03 which cor- 
responds to Z = 0.31 + 0.04 Z . We observed no significant gradients in 7\Oiii), n e or [O/H] 
on a scale of 73 pc extending in four directions from the ionising source. The properties of 
the H ii region can be reproduced by a CLOUDY model which uses the central cluster as ion- 
ising source, thus providing a self-consistent interpretation of the data. We also report on the 
serendipitous discovery of He n nebular emission associated with the nearby source NGC 55 
C2_35, a feature usually associated with strong X-ray sources. 

Key words: stars: early type - stars: massive - stars: Wolf-Rayet - stars: individual: NGC 55 
Cl_31 - galaxies: ISM - galaxies: individual: NGC 55 



1 INTRODUCTION 

The most luminous stars in low metallicity galaxies are of special 
interest. These may have masses exceeding 200 M Q , defining the 
upper mass limit of stars. Until recently, such massi ve objects were 
only found in the cores of young massive clusters JCrowther et alj 
120101) . but the first such object has now been found in apparent 
isolation dBestenlehner et al] l201ll) . This motivates a search 
for very massive stars in Local Group dwarf galaxies, or even 
more distant systems, where current instrumentation does not 
yet allow stellar clusters to be spatially resolved. Depending on 
mass and metallicity, the mass-loss rates of the brightest stars 
may be so high that their winds become optically thick, resultin g 
in hydrogen-rich WN spectra dde Koter, Heap & Hubenvl 1 19971) . 
Such t argets provide impor t ant tests for the theory of line-driven 
winds (Grafener et al. 20lil lvirik et a i] |201lh . Massive stars up to 
about 80 M in metal poor environments receive special attention 



as well, as those that have a rapidly rotating core at the end of their 
lives may produce broad line Type Ic supernovae or hypernovae, 
which are perhaps connected to long duration gamma-ray bursts 
dMoriva et aifcoid) . 

In this context, massive stars in the Magellanic Cloud s 
have been extensively studied (see e.g. lEvans et alj|2004l 1201 lh . 
Extending such studies to more distant galaxies requires sensitive 
spectrographs mounted at the largest telescopes and, so far, has 
only be attempted at low resolution l lBresolin et alj I2006L 120071 : 
ICastro et alj|2008l) . Such low resolution studies can be hampered 
by the presence of nebular emission, as due to their strong UV 
flux, massive stars ionise the ambient environment creating Hn 
regions. Through observations at higher spectral resolution, 
nebular emission, rather than being a complicating factor, may 
help to further constrain physical pr operties of the ionising source 
(see e.g. lKudritzki & Hummer! 199Ch . 



* Based on observations made with ESO Telescopes at the Paranal Obser- 
vatory under program 60A-94 19(A) 
t O.E.Hartoog@uva.nl 



In this paper we take the first step towards quantitative spec- 
troscopy of massive stars outside the Local Group. We present the 
first medium-resolution spectrum of a luminous early-type source 
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Figure l.B-V-Ha image of the host galaxy NGC 55 (MPG/ESO, Wide 
Field Imager) with a zoomed in region (FORS Ha image) that indicates the 
location of our source. The rectangles show the projections of the VIS en- 
trance slit (0.9" X 11") at the different observing epochs: (1) is 13/08/09, 
(2) is 27/09/09 and (3) is 30/09/09. C1.31 is located where the slit projec- 
tions cross. C2_35 is another massive star candidate, of which we see the 
surrounding nebular emission in our slit (see Section l4~2V 



in NGC 55 (~2.0Mpc) and its surrounding region, obtained with 
the new X-shooter spectrograph at the ESO Very Large Telescope 
(VLT). The medium spectral resolution (R ~ 6000) and unique 
spectral coverage of X-shooter allow for a detailed analysis of the 
stellar spectrum, and, importantly, for an improved correction of 
the nebular emission lines that can be distinguished from the un- 
derlying stellar spectrum. 



NGC 55 

NGC 55 (see Fig. [TJ is located in the foreground of the Sculptor 
Group , at a distance of approximately 2.0 Mpc dGieren et alj 
I2008L and refererences therein). Though it is difficult to deter- 
mine its morphological type due to i ts high inclination of ~80° 
jHummel. Dettmar &Wielebinskl [l9861. the galaxy is likely of 
type SB(s)m: a barred spiral with an irregular appearance, very 
simil ar to the Large Magellanic Cloud dde Vaucouleurs & Freeman] 
Metallicity measurements of NGC 55 show a range of 
values between 0.23 and 0.7 Z Q , all det ermined by analy sis of 
forbidden oxyg en line emiss ion (se e e.g. IWebster & Smithll 19831; 
Stasihska, Comte & Vigroux) 1 19861 : IZaritskv. Kennicutt & Huchral 
19941 : iTullmann et alJl2003jT 

The blue massive star population of NGC 55 has bee n studied 
in the context o f the Araucaria proj ect jGieren et al.ll2005b . As part 
of this project. ICastro et alj ( 2008) have presented low resolution 
(R = 780) optical (390 - 490 nm) VLT/FORS2 spectra of approx- 
imately 200 blue massive stars in NGC 55. In search for the most 
massive star in this galaxy, we selected NGC55C1_31 (for the 
remainder of this paper Cl_31, RA 0:15:00.01, DEC -39:12:41.39, 



Table 1. Overview of the observations. Date and mid-exposure time of the 
observations; Seeing measured from R band acquisition images; Angular 
distance from the source to the Moon and fraction of lunar illumination 
(FLI); Exposure time per arm; Position angle (North to East) of the slit on 
the sky (along the parallactic angle at the time of the observations). 





Date 


Seeing 


Moon dist. 


Exp. T. 


Pos. Ang. 




Time (UT) 


(R-band) 


FLI 


(s) 




1 


13/08/09 


0.69 - 0.88" 


69° 


2x900 


53.4° 




08:52 




56% 






2 


27/09/09 


0.87 - 0.93" 


66° 


4x900 


-52.3° 




03:43 




60% 






3 


30/09/09 


0.96- 1.50" 


45° 


4x900 


-56.4° 




03:26 




85% 







Table 2. Overview of the X-shooter instrument properties. Instrument arm; 
Wavelength range; Projected slit size; Measured resolving power R = A./AA; 
Resolving power according to the X-Shooter User Manual. 



Arm 


Range (nm) 


Slit dimensions 


R 


Rth 


UVB 


300 - 560 


0.8" X 11" 


6268 ± 179 


6200 


VIS 


550 - 1020 


0.9" x 11" 


7778 ± 264 


8800 


NIR 


1020 - 2480 


0.6" x 11" 


7650 ± 258 


8100 



indicated in Fig. Q3 because of its brightness and its classification 
as an early O-type supergiant. 

In the following section we describe the observations and data 
reduction for both the stellar spectrum Q2.lt and the nebular spec- 
tra as a function of location along the slit Q2.2t , In Section|3] we 
constrain the overall properties of the central source by compar- 
ing the observed hydrogen and helium line profiles of Cl_31 with 
simulated profiles. In Section [4] we analyse the nebular spectra. In 
Section [5] we explore the effect of luminosity and temperature of 
the ionizing source on the properties of the surrounding nebula, re- 
sulting in a consistent picture on both the properties of the central 
ionizing source and the surrounding nebula (^6}. We summarise the 
main conclusions in Section|7] 



2 OBSERVATIONS AND DATA REDUCTION 

The observations of Cl_31 were obtained as part of X-shooter 
Science Verification (SV) Runs 1 and 2. Spread over three nights, 
the total exposure time is 2.5h. The observations were carried out 
in nodding mode using a nod throw of 5". We refer to Table[T]for 
the details of the observations and the observing conditions. 

Light that enters X-shooter is split in three arms using 
dichroics: UV-Blue (UVB), VISual (VIS) and Near-IR (NIR). Each 
instru ment arm is a fixed format cross-disperse d echelle spectro- 
graph dD'Odorico et alj2006l : IVernet et alj201 ll) . Table|2]gives for 
each arm the wavelength range, the projected dimensions of the slit 
and the resolving pow er. The data are reduced with the X-sh ooter 
pipeline version 1.2.2 ( Modigliani et ai] |2010l : [Goidonill201lh . Al- 
though our source is observed in nodding mode, we have reduced 
the UVB and VIS science frames separately, using the staring mode 
reduction recipe. We follow the full cascade of X-shooter pipeline 
steps ('physical model mode'), up to obtaining two-dimensional 
(2D) straightened spectra, without sky subtraction. See Section ETI 
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Table 3. Signal-to-noise ratio (SNR) per resolution element of different 
parts of the UVB and VIS spectra per night, and of the final combined 
spectra. 



UVB arm 
Night 


[424:428] 


SNR [range (nm)] 

[460:465] [505:510] 




13/08/09 
27/09/09 
30/09/09 


16.9 
14.0 
7.2 


21.3 
21.2 
9.8 




21.4 
13.7 
9.8 




combined 


21.2 


29.3 




27.4 




VIS arm 
Night 


[604:609] 


SNR [range (nm)] 
[675:680] [811 


816] 


[975:978] 


13/08/09 
27/09/09 
30/09/09 


12.1 
10.3 
6.2 


15.5 
12.0 
9.8 




22.0 
19.8 
13.3 


6.3 
5.2 
4.8 


combined 


16.2 


20.9 




30.0 


8.7 



for the steps through which we obtain the stellar spectrum of C 1 _3 1 , 
and Section |2~2l for the extraction procedure of the nebular spectra 
along the spatial direction of the slit. 



2.1 The sky-corrected stellar spectrum 

UVB and VIS arms 

One-dimensional object and sky spectra are extracted from the 2D 
spectra. The sky spectrum is extracted from regions with the low- 
est continuum and nebular line emission. In the nights with four 
consecutive exposures, cosmic rays have been removed by taking 
the median value of each pixel in the four exposures. The sky spec- 
tra are subtracted from the object spectra, thereby correcting for the 
contamination by moonlight as well. After applying the barycentric 
correction, the sky-corrected object spectra of the three nights have 
been combined. Table [3] gives the signal-to-noise ratio (SNR) per 
resolution element of the result of each night and of the combined 
spectrum. For the stellar line analysis we normalised the spectrum. 
Independently, we also calibrated the flux of the UVB and VIS 
spectra with standard star BD+17° 4708 (sdF8) taken on 27/09/09, 
since there were no appropriate flux standard observations for each 
individual night. We did not correct the spectrum for slit losses. 



NIR arm 

We reduced the NIR spectra following the pipeline cascade for nod- 
ding mode up to extracted ID spectra, and we combined the three 
nights. A telluric standard star (HD4670, B9 V) is used to cor- 
rect for the telluric absorption features in the combined NIR spec- 
trum, and to calibrate the relative flux. The NIR spectrum is scaled 
to match the absolutely calibrated VIS spectrum. The SNR of the 
combined NIR spectrum of all three nights is ~ 3 in the J band, and 
< 1 in the H and K bands. This is too low to detect stellar or neb- 
ular features in the spectrum; however, the level of the continuum 
can be retrieved by binning the flux in the atmospheric bands. The 
result is shown in Fig. [2] 



2.2 The nebular emission spectra 

The 2D spectra of each night are combined using the median 
value. We flux calibrate the combined 2D spectra with the same 
photometric standard star as we used for the object spectrum. Fol- 
lowing the trace of Cl_31, we extract 20 sub-apertures from each 
night's combined 2D spectrum, resulting in a set of ID nebular 
spectra that are spatially separated by 0.75" eaclQ. Independently, 
we apply this combining and sub-aperture extraction procedure 
as well to the non-flat-fielded 2D spectra, now using the sum. 
This allows us to determine the number of photons N, and thus 
the photon noise error for every emission line in the nebular 
spectra. These errors are propagated in the values of the nebular 
properties described in Section|4] 



3 ANALYSIS: THE C1.31 STELLAR SPECTRUM 

Fig. shows the combined flux-calibrated stellar spectrum of 
Cl_31. The extinction corrected flux should show a Rayleigh- 
Jeans wavelength dependence (F A cc A~ 4 ), b ecause we expect a n 
early-type star based on the classification in ICastro et al.1 (2008). 
To lift the observed spectrum to the slope of the scaled Kurucz 
model of a B0 star, we need to de- redden our spectr um with 
Av, 8tar = 2.3 ± 0.1, adopting R v = 3.24 (Gieren et al. 2008). We ap- 
ply the parametrized extinction law by Cardelli, Clayton & Mathisi 
(1989). As a check, we do the same exercise by varying both R v 
and A Vfitm . With 3.0 < i?v < 3.5, the spectrum can be de-reddened 
to the intrinsic F A cc A slope with an Ay; star = 2.3 ±0.1. With R v 
outside this range, the de-reddened spectrum does not match the 
models for any value of Ay star . 



ICastro et ai] d2008l) report magnitudes V = 18.523 and 
/ = 19.239 for NGC55C131, which were obtained as part o f 
the Araucaria Cepheid search project (Pietrzvhski et al.1 hood) . 
The / magnitude is in good agreement with our flux-calibrated 
spectrum; V is not. The reported color V - I = -0.716 is even 
bluer than a theoretical Rayleigh- Jeans tail, suggesting a problem 
with the photometrjQ. From our flux-calibrated spectrum we 
obtain V =* 20.1 ± 0.1. Using A Kstar = 2.3, and d = 2.0 Mpc, 
we derive an absolute magnitude of My = -8.7±0.4 for the source. 

Fig. [3] compares the observed UVB spectrum of Cl_31 
between 380 and 490 nm w ith that of spectral standard stars 
dWalborn & Fitzpatrickiri99rj|) . In the Cl_31 spectrum the B aimer 
and He i lines show artifacts of the nebular emission correction. 
However, most of the line wings are left unaffected thanks to 
the relatively high spectral resolution of X-shooter. Based on the 
non-detection of Hen/14541 and the presence of He 1/14471 in the 
Cl_31 spectrum one can not classify this source as an early O-star. 
By adding artificial noise to the standard star spectra, matching 
the SNR of our observations in this range, we estimate that the 
He H/14541 line of a star with spectral type later than 07.5 will 
not be detectable, thus suggesting a late spectral subtype. In Sec- 
tion [37TTT] we will use stellar atmosphere models to constrain the 



1 The sub-aperture size of 0.75" is just below the average FWHM of the 
seeing. 

2 This has been confirmed by the authors. Corrected photometric values are 
V = 19.87 ± 0.05, / = 19.25 ± 0.05 (priv. comm. Pietrzyfiski), in agreement 
with our findings. 
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Figure 2. Flux-calibrated spectrum of NGC55 C 1_3 1 (bottom), and the same spectrum (top) corrected for extinction with Ry = 3.24 and A^star = 2.3 (see 
text). The dashed li ne shows the FAS TWIND spectral energy distribution for the best fitting combination of models (see Section|6j; the solid line is a Kurucz 
model of a B0 star <Kuruczll979lll993l) . The ranges of X-shooters instrument arms are indicated at the top of the graph. The open squares point to the updated 
V and / magnitudes (prfv. comm. Pietrzyhski). For clarity, the NIR spectrum is shown here with a smoothing of 7 points. The triangles show the integrated 
values of the observed flux in the J and H bands; the diamonds show their extinction corrected values. 
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Figure 3. The normalised sky-corrected object spectrum (top) between 380 
and 490 nm compa red with the spec tra of standard O supergiants from 03 to 
Q8.5 iWalborn & Fitzpatrick 1990). Note the very broad Heii/14686 emis- 



effective temperature more quantitatively. Supergiants may have 
strong emission in Heii/14686 due to their stellar wind, but none 
of the standard stars show a feature as broad as that in our observa- 
tions. Broad Heii/14686 emission lines are the strongest features in 
Wolf-Rayet stars of type WN. This line profile will be analysed in 
more detail in Section [3.1.4l 

When we compa re the X-shooter sp ectrum to the FORS2 
spectrum of Cl_31 in ICastro et al] |2008), we conclude that the 
general appe arance is similar, including the shape of He n /14686. 
ICastro et al] observed Heii/i4200 weakly in emission as well; this 
we can not confirm. The FORS2 observations will have suffered 
from similar nebular contamination in the Balmer lines as well as 
in some He i lines, which may have hampered earlier classification, 
but correcting for this is even more difficult at lower spectral reso- 
lution. 



3.1 Stellar line profile modelling 

Although we are going to make clear in this paper that Cl_31 is 
very likely a composite source, we first approach the spectrum as 
if it is produced by a single star. The derived physical parameters 
therefore represent averages of the flux-weighted components 
contributing to the spectrum. We note, though, that the integrated 
light from clusters - in particular the hydrogen ionising radiation 
- is often dominated by only a few of the most massive components. 

The profiles of spectral lines are responsive to various stellar 
parameters such as effective temperature, mass-loss rate, surface 
gravity, che mical abundances and rotation speed. We have used 
FASTWIND JPuls et al]|2005h to model stellar atmospheres and to 
compute profiles of spectral lines. FASTWIND calculates non-LTE 
line blanketed stellar atmospheres and is suited to model stars with 
strong winds. We first tried to apply a genetic fitting algorithm 
with FASTWIND models (see lMokiem et alj20oA to the observed 
spectrum, in order to fit a large number of parameters at the same 
time. This did not result in well constrained parameters, because 
the shape and width of He n /i4686 could not be fitted at the same 
time as the other lines. 

Instead we constructed a grid of FASTWIND models, and 
constrained the parameters by comparing the observed profiles 
with the models. The main grid varies effective temperature 
T eff = 27 500 - 35 000 K with steps of 2 500 K; and covers values 
for the mass-loss rate M = 1,3,6 and 10 x 10~ 6 M Q yr~' and 
luminosity log(L,/L G ) = 5.4,5.6 and 5.8. All models have a 
mass M = 40 M Q , wind acceleration parameter /3 = 1.0, and 
helium to hydrogen number density Nn c /Nn = 0.1. Radius R, 
and surface gravity g are fixed by the other parameters. The 
terminal wind velocit y is assumed to be 2.6 times the surface 
escape velocity (see lLamers. Snow & Lindholm|[l995r) . In Table g] 
we show the parameter values of a subset of the grid, i.e. the 
models discussed in more detail in this paper. Unless stated 
otherwise, the synthesized line profiles are produced using a 
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Table 4. Parameters used for the described FASTWIND models, which are a subset of a larger grid. R t , log g and v«, = 2.6 X v esc depend on the other variables. 
The first four models have the same M, L, and M, but different T e g causing the remaining parameters to be modified. The next three models are variations on 
MOD3 1 , but with different values for M. MOD69 is like MOD3 1 , but with an artificially high wind acceleration parameter ft. The last two models MOD89 
and MOD90 provide the best fitting multiple model configuration (in a visual flux ratio ~ 4 : 1). The bolometric luminosities L* are chosen such that they 
reproduce the measured My with ten times MOD89 and two times MOD90. 



ID 


M 


T e ff 


L, 




R, 


M 


Voo 






z 




(Mo) 


(K) 


(L ) 


(cm s~ 2 ) 


(R e ) 


(M yr-') 


(kms^ 1 ) 






(Z Q ) 


MOD39 


40 


27 500 


io 5 - 6 


2.99 


33.69 


6.00 x 10~ 6 


1750.33 


1.0 


0.1 


0.3 


MOD31 


40 


30000 


10 5.6 


3.30 


23.40 


6.00 x 10~ 6 


2100.40 


1.0 


0.1 


0.3 


MOD32 


40 


32 500 


10 5.6 


3.44 


29.94 


6.00 x 10~ 6 


2275.43 


1.0 


0.1 


0.3 


MOD33 


40 


35 000 


io 5 - 6 


3.57 


17.19 


6.00 x 10~ 6 


2450.47 


1.0 


0.1 


0.3 


MOD10 


40 


30000 


io 5 - 6 


3.30 


23.40 


1.00 x 10~ 6 


2100.40 


1.0 


0.1 


0.3 


MOD11 


40 


30000 


io 5 - 6 


3.30 


23.40 


3.00 x IO" 6 


2100.40 


1.0 


0.1 


0.3 


MOD 12 


40 


30000 


io 5 - 6 


3.30 


23.40 


1.00 x 10~ 5 


2100.40 


1.0 


0.1 


0.3 


MOD69 


40 


30000 


io 5 - 6 


3.30 


23.40 


6.00 x IO" 6 


2100.40 


3.0 


0.1 


0.3 


MOD89 


30 


30000 


10 5.24 


3.54 


15.46 


2.00 x 10~ 6 


2237.86 


1.0 


0.1 


0.3 


MOD90 


80 


50000 


10 583 


4.26 


10.98 


3.00 x 10~ 5 


4336.76 


1.0 


0.8 


0.3 
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Figure 4. Observed line profiles compared with model line profiles com- 
puted for various effective surface temperatures (see text). 



files of Heii4387, Hei^4471, Hen^4200 and HeiU4541 along 
with profiles from atmosphere models that only differ in effective 
temperature. A visual comparison of the profiles shows that the 
models with r cff < 35 000 K best reproduce the spectrum as oth- 
erwise Hen would have been detected, and He 1/14387 would not 
have been as deep. 



3.1.2 Rotational velocity 

Since the observed profile of He i /14387 is not affected by nebu- 
lar emission, we can use it to constrain the characteristic rotational 
broadening. Fig.|5]shows the HeU4387 profile of MOD31 (T eS = 
30 000 K) convolved with rotational profiles to simulate three dif- 
ferent rotational velocities: v ro tsin(i) = 50, 150 and 250 km s" 1 . 
The line with v rot sin(i') = 50 km s _I is clearly not broad enough to 
fit the observed profile, while the line with v rot sin(i) = 250 km s -1 
appears to be too broad. We conclude that the width of the lines is 
best reproduced by models with v rot sin(i) = 150 + 50 km s . 



microturbulent velocity v lur b = 10kms~' and projected rotational 
velocity v lot sin(f') = 150 km s _I , with i the inclination of the stellar 
rotation axis relative to the plane of the sky. An instrumental 
profile matching the resolution of X-shooter is applied as well. 

We first investigate the overall impact of the effective tem- 
perature Q3.1.U . rotational velocity Q3.1.2t . and mass-loss rate 
(j 33.1.3t by comparing with models involving a single star. In 93.1.41 
we show that the observed profile of He n /14686 can not be repro- 
duced with a single star. In §|3.1.5| we will discuss the luminosity 
constraint that follows from My. 

3.1.1 Effective temperature 

The He i and He n lines can be used to determine the characteristic 
effective temperature of the source. Fig.|4]shows the observed pro- 



3.1.3 Mass-loss rate 

In the nebular spectrum He is very strongly in emission. This is 
mostly due to the surrounding Hn region. In the sky-corrected 
object spectrum (Fig. [6}, the Ha- line wings are clearly visible, and 
reveal that Ho- is in emission in the spectrum of Cl_31 as well. 

The profile of Ha is very sensitive to the mass-loss flux 
M/4nRl. In the right panel in Fig. [6] we plot the observed spec- 
trum together with the line profiles from single-star models with 
r e ff = 30000 K and R, = 23.4 Rq, and four different mass-loss 
rates: 1 x 10~ 5 , 6 x 10~ 6 , 3 x 10~ 6 , and 1 x 10~ 6 M Q yr 1 . De- 
pending on the mass-loss rate, the line is either in emission or 
in absorption. The computed profile for MOD31 reproduces the 
observed profile best, and corresponds to a mass-loss flux close 
to ~9 x 10~ 10 M yr 1 R- 2 (i. e. M = 6 x 10~ 6 M Q yr 1 for 
R, = 23.4 Rq). 
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Figure 5. The observed line profile of He i /14387 compared with the same 
model profile computed for various projected rotational velocities. 
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Figure 6. Observed line profiles compared with model line profiles for var- 
ious mass-loss rates (left panel) and projected rotational velocities (right 
panel). 



3.1.4 He ii A4686 

In the spectrum of Cl_31, Hen/i4686 is in emission, which is 
a common feature in spectra of O-type supergiants (see Fig. [3}- 
However, in none of the comparison spectra this line is as broad 
as in our spectrum (~3000 km s~'). The equivalent width is only 
-3.6 ± 0.04 A, i.e. much weaker than in typical WN star spectra. 
The feature is of stellar origin as no nebular counterpart is detected. 

The strength of He n /14686 depends on various parameters in 
the FASTWIND models. The modeled line can be made stronger by 
(1) increasing the effective temperature (a larger fraction of the he- 
lium will be ionised), (2) increasing the mass-loss rate, (3) increas- 
ing the helium abundance, (4) increasing the wind acceleration pa- 
rameter p or (4) decreasing the terminal wind velocity. But these 
modifications only make the line stronger, not much broader. To 
simulate the profile of the line in both strength and width using rota- 
tional broadening, a rotational velocity of v rot sin(i) = 1200 km s 
is needed: Fig. [6] shows one of the models (MOD69) for which 
this line has an equivalent width of -3.58 A, convolved with a ro- 
tational profile simulating rotational velocities of v rot sin(;) = 150, 
600, 900 and 1200 kms -1 . Only when v rot sin(0 ~ 1200 km s -1 , the 
line is broad enough to reproduce the observed profile. This v rot is, 
for acceptable radii and masses, higher than the escape velocity. 
He ii /14686 is formed in the wind, therefore its width can be much 
higher than the surface rotational velocity. But with such a dense 
and fast wind, Ho- would have been much broader as well, and the 
photospheric He i and He n absorption lines would not be visible. 
Therefore we conclude that the Heii/14686 emission line has a dif- 
ferent origin: it is a diluted WN feature (see Section|6]l. 



Table 5. Reddening corrected (see text) line flux ratios with respect to H/j 
for the strongest unblended emission lines in the nebular spectrum close to 
our central source. 



Line 








13/08/09 
Ratio 


27/09/09 
Ratio 


Average 
Ratio 


[On] 


A 


3726 





1.554 ±0.050 


1.476 ±0.033 


1.500 ±0.028 


[On] 


A 


3728 


8 


2.271 ±0.065 


2.188 ±0.044 


2.214 ±0.036 


H-9 


A 


3835 


4 


0.062 ±0.011 


0.060 ± 0.007 


0.061 ±0.006 


[Nenfji 


3868 


8 


0.325 ±0.018 


0.306 ±0.012 


0.312±0.010 


H-8 


A 


3889 


1 


0.201 ±0.014 


0.198 ±0.009 


0.199 ±0.008 


Hd 


A 


4102 





0.275 ±0.014 


0.271 ±0.010 


0.272 ±0.008 


Hy 


A 


4340 


5 


0.517 ±0.019 


0.505 ±0.013 


0.509 ±0.011 


[O hi] 


A 


4363 


2 


0.039 ± 0.004 


0.037 ± 0.003 


0.037 ± 0.002 


He i 


A 


4471 





0.034 ±0.005 


0.034 ± 0.004 


0.034 ± 0.003 


Hen 


A 


4686 





0.006 ± 0.002 


0.003 ± 0.001 


0.003 ± 0.001 


H/j 


A 


4861 


3 


1.000 ±0.027 


1.000 ±0.019 


1.000 ±0.016 


[Out] 


A 


4958 


9 


1.171 ±0.029 


1.144 ±0.020 


1.153 ± 0.016 


[Out] 


A 


5006 


7 


3.477 ±0.074 


3.426 ±0.052 


3.443 ±0.043 


He i 


A 


5876 





0.1 10 ±0.007 


0.1 13 ±0.005 


0.1 12 ±0.004 


[Sm] 


A 


6312 





0.015 ±0.002 


0.017 ±0.001 


0.016 ±0.001 


[Nil] 


A 


6548 





U.UOU ± U.UU4 


o _l n om 


U.Uol ± U.UU2 


Ha 


A 


6562 


8 


3.040 ±0.063 


3.052 ± 0.045 


3.048 ±0.037 


[Nn] 


A 


6583 


4 


0.179 ±0.007 


0.193 ±0.005 


0.188 ±0.004 


He i 


A 


6678 





0.029 ± 0.003 


0.028 ± 0.002 


0.028 ± 0.002 


[Sn] 


A 


6716 


5 


0.269 ± 0.009 


0.292 ± 0.007 


0.283 ± 0.005 


[Sn] 


A 


6730 


8 


0.189 ±0.007 


0.206 ± 0.005 


0.200 ±0.004 


[Arv] 


A 


7005 


9 


0.009 ± 0.001 


0.006 ± 0.001 


0.007 ±0.001 


[Arm] 


A 


7135 


8 


0.087 ± 0.004 


0.086 ± 0.003 


0.087 ±0.002 


[Arm; 


A 


7751 


1 


0.022 ± 0.001 


0.022 ± 0.001 


0.022 ± 0.001 


Pa- 10 


A 


9015 





0.017 ±0.001 


0.018 ±0.001 


0.018 ±0.001 


[Sm] 


A 


9068 


9 


0.200 ±0.005 


0.203 ± 0.004 


0.202 ±0.003 


Pa-9 


A 


9229 





0.025 ± 0.002 


0.024 ± 0.001 


0.024 ±0.001 


[Sm] 


A 


9531 





0.475 ±0.011 


0.437 ± 0.008 


0.449 ± 0.006 


Pa-7 


,110049 


4 


0.045 ± 0.004 


0.046 ± 0.003 


0.046 ± 0.002 



3.1.5 Luminosity 

In the calibration of O-stars bv lMartins. Schaerer & Hillierl J2005h 
the visually brightest supergiant has M v = -6.35. Cl_31 is with 
M v = -8.7 almost an order of magnitude brighter. It is therefore 
likely that Cl_31 is a composite object such as a cluster containing 
several luminous stars. The slit width (0.8" for UVB) at a distance 
of 2.0 Mpc corresponds to a physical size of 7.8 pc, which is large 
enough to contain for example the Orion Trap ezium cluster, or 
even more ma ssive open cluste rs such as Tr 14 dSana et alj[201oh 
or NGC 6231 dSana et alj|2008l) . 

In summary, the normalised line profiles of hydrogen and he- 
lium, except He n /14686, can be reproduced by a single-star model 
with parameters 7/ eft ~ 30000 K and M ~ 6 x 10~ 6 M Q yr -1 , 
R, ~ 23.4 R Q and v rot sin(0 = 150 ± 50 km s -1 , i.e. a late O su- 
pergiant star (MOD31 in Table [4j. A star with this temperature 
would need a luminosity of log(L„/L G ) ~ 6.3 to reproduce M v , 
which is too high for a single O-type supergiant. The considera- 
tions regarding the Heii/14686 line also point in the direction of 
the spectrum being a composite of different sources, weighed by 
their visual brightness. This scenario will be explored in Section[6] 



4 ANALYSIS: THE NEBULAR EMISSION LINE 
SPECTRUM 

The nebular emission spectra show hydrogen recombination lines 
of the Balmer and Paschen series, Hei lines, and forbidden lines 
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of On, Om, Sn, Sin, Nn, Nem, Arm, Ariv and Arv. To these 
spectra, no sky-correction could be applied, because nebular 
emission covers the full slit. Therefore, every nebular line we 
measure might have a contribution from the sky continuum. 
We minimize this error by subtracting the local continuum next 
to the line in wavelength. The ratios with respect to H/3 of the 
nebular lines at the position of our source are listed in Table [5] 
The extinction corrected integrated specific intensity of H/3 is 
5.8 x 10~ 15 erg s -1 cm" 2 arcsec~ 2 . 

We use the diagnostics for e lectron temperature and oxygen 
abundance from lPagel et"all dl992l) : 



T = 



1.432 



log/? - 0.85 + 0.03 log T + log (1 + 0.0433xr 006 ) 



(1) 



where T = T (O m), the electron temperature in the region where 
oxygen is doubly ionized, in units of 10 4 K. R and x are given by 



Al4959 + ^J5 



A)4363 

x = 10~ 4 ;i c 77 1/2 , 



(2) 
(3) 



where / is the integrated specific intensity intensity of the indicated 
emission line and n e is the electron density in cmT 3 . Ti is the elec- 
tron temperature in units of 10 4 K i n the singly ioniz ed region, and 
follows from model calculations bv lStasinskal dl990h : 

r 2 - ] = [r(On,Nii,Sii)r l = 0.5(7-' +0.8) . (4) 

The mean ionic abundance ratios for O n and O m along the line of 
sight are calculated as follows: 

,„ , { f~\ / t t \ 1 ^.13726 +^.13729 - „„„ 1.676 

12 + log(On/Hn) = log +5.890+ 

-0.40 log T 2 + log (1 + 1.35x); (5) 

1/1 i //-v /tt \ i Al4959 + ^15007 , 1.251 

12 + log(Om/Hii) = log : +6.174 + 



-0.55 log 7/; 



(6) 



The total oxygen abundance is obtained by adding Eqs. [5] and [6] 
i.e. by assuming these two ionisation stages to be dominant in the 
H ii region. 



4.1 Nebular emission properties along the slit 

Fig. [7] shows the nebular properties along the slit for the spectra 
obtained on 13/08/09 and 27/09/09. The orientation of the slit dif- 
fers between these observations: positive offset is ~North-East for 
13/08/09 and -North- West for 27/09/09 (see Fig. QJ. We do not 
show the results for the 30/09/09 observation, which has a similar 
position angle and gives a similar result to 27/09/09, though with 
larger error bars. As mentioned in Section |2"7I1 the errors are ob- 
tained by propagating the photon noise on the intensity of the used 
lines. We choose n c = 20 cmr 3 arbitrarily, but in agreement with the 
low density limit of the density-sensitive [O n] and [S n] ratios we 
measure (see panel (d) in Fig.[7j and the analysis below). The error 
on the density is not propagated into the errors on the other parame- 
ters, because they all depend very weakly on the density (see Eqs.Q] 
and[3}. In the following sections we discuss the panels in Fig.[7j 



Stellar continuum and nebular emission 

Panel (a) shows the shape of the continuum between 411.0 - 
433.5 nm; here we see the trace of our central source. We do not 
take into account the difference in A v along the slit (see panel c). 
The nebular emission in H/3 and [O m] /15007, two of the strongest 
emission lines in the spectra, is also shown. [O m] /15007 is a for- 
bidden transition and thus only emitted by low density nebulae. The 
27/09/09 observation suggests that the nebular emission is built up 
from various discrete peaks, of which one is centered on Cl_31. 
This part of the nebula, with a radius of ~20-30 pc, is likely ionized 
by Cl_31. Therefore, the nebular properties at the central location 
in the slit are related to the properties of Cl_31 (See Section|5]l- 



HeuA4686/Hf} 

Panel (b) shows the ratio of the He n /14686 and H/3 nebular emis- 
sion line. Around Cl_31, this nebular line is not present, but it is 
very pronounced in the 13/08/09 spectrum around offset -4.5". 
This is also clearly visible in the 2D spectrum around this line (top 
left image of Fig. |7j. This feature will be discussed in more de- 
tail in Section l4~2l This nebular feature is much narrower than the 
wind He H/14686 line we discussed in Section [3. 1.41 The nebular 
He n /14686 line has a Gaussian FWHM of ~ 1 A like the other neb- 
ular lines, slightly larger than the resolution element at this wave- 
length (-0.75 A). 



Extinction 

Panel (c) shows the values of A Vncb , which is the extinction de- 
rived from the IS M hydrogen line ratios. Adopting R v = 3.24 
dGierenetalj2008l) .A v,neb is found for each aperture by minimizing 
the following expression: 



2 



v v i-t 



(7) 



with qo the theoretical ratio for Case B rec ombination in the low 
densi ty limit at T = 10 000 K (see e.g. lOsterbrock & Ferland 
2006), q the inte nsity of a line with respect to H/3 after applying 
A v ,neb (following ICardelli. Clayton & Mafhis! 1 19891) . cr 2 the vari- 
ance on q, and v the number of degrees of freedom. We use H-9, 
Hy, H<5, Ho-, Pa- 10, Pa-9 and Pa-7. The confidence interval on 
Ay.neb is given by the value for which x 2 mA rises by 1- ^v.ncb is used 
in the derivation of the properties per aperture shown in panels (d) 
to (f) of Fig. UJ the error on A VjK b is not propagated as its influence 
on the other parameters is small. 

Av, n <;b at the location of Cl_31 is derived to be 1.30 + 0.15 and 
1.15 ± 0.10 in the 13/08/09 and 27/09/09 spectra respectively. 
This is lower than A VtStm = 2.3 ± 0.1 (see Section [3] and Fig. [2}. 
Furthermore, w e note that both Ay ne b and Ay, st ar are larger than 
A v = 0.45 from lGieren et alj d2008l) . However, this latter value is 
an average for NGC 55 as a whole. This range of values reflects 
local variations and are to be expected, especially in an almost 
edge-on galaxy. 



Electron density 

Panel (d) shows [O n] and [S n] ratios that are sensitive to electron 
density (see e.g. lOsterbrock & Ferland 120061) . For a value > 1.4, 
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Figure 7. Spatial profiles of emission lines and properties measured from the emission line spectrum in 20 sub-apertures along the slit On the x-axis is distance 
in ai'csec from Cl_31. The left and right frames, labeled by their observing dates, correspond to observations with different slit orientations, see Fig.[T]and 
Table[T] The image at the top is a part of the 2D spectrum around 468. 6nm. Panel (a) shows the average level of the continuum between 411.0 - 433.5 nm, 
without taking into account the difference in Ay nc b along the slit. It also shows the intensity of the two strong nebular lines H/J and [O in] /15007, with respect to 
the value at the location of the central source. Panel (b) shows the flux of the He n /t4686 nebular emission line with respect to H/3. Panel (c) shows the value of 
Ay ne jj, derived from a number of Balmer and Paschen emission lines, assuming Case B recombination. Panel (d) shows [O ii] and [S n] ratios that are sensitive 
to electron density. Panel (e) shows T (O m) calculated from the ratio R of the [O m] lines (Eq.[T}. Panel (f) shows the total oxygen abundance 12 + log (O/H). 



both ratios are in the low density (< 10 2 cirr 3 ) limit, so this mea- 
surement only provides an upper limit. Because these line pairs are 
very close in wavelength, their ratios are not affected by Ay iIle b. 

Electron temperature 

Panel (e) gives the electron temperature calculated from the ratio R 
of the [O m] lines (EquationQ~|l. In the 27/09/09 observation (Fig. [7] 
right) in the apertures with offset > 2", [Om] /14363 was hardly 
detected resulting in an underestimation of the error in T (O m). 
Taking this into account, and weighing the better quality spec- 
tra more strongly, we conclude that around the location of Cl_31 
T (O m) = 1 1500 ± 600 K, and that there are no significant gradi- 
ents in the two spatial directions indicated in Fig.Q] 



Oxygen abundance 

Panel (f) shows the total oxygen abundance [O/H] = 12 + 
log (O/H). The slight drop in [O/H] that we see at offset > 
2" in the 27/09/09 observation is a propagated effect from 
the uncertain determination of r(Ora). Excluding this region, 



we find an average of [O/H] = 8.18 ± 0.03, which corre- 
sponds to_Z_^a31± 0.04 Z Q adopting [0/H] = 8.69 ± 
0.05 dAsplund et all [2009). Though low er oxygen abundance s 
are reported for NGC55 (8.08 ± 0.10, iTullmann etal] l2003h . 
on average slightly hi gher v alues are measured (8.23 - 8.39 , 
IWebster & Smith! 1 19831: 8.53. IStasinska. Comte & Vigrouxlll986l 



8.35 ± 0.07. IZaritskv, Kennicutt & Huchrl l 994 ) . 



4.2 C2_35 

C2_35 (RA 0:14:59.68, DEC -39:12:42.84) is located 4.5" West 
South-West of Cl_31 (see Fig. [I). Its strongly ionised surrounding 
nebula is visible in the 13/08/09 spectrum (Fig. [7] left, panel a). 
At offset -4.5" we see a weak continuum, as the point-spread 
function of C2_35 is mostly outside the slit. We detect strong 
(forbidden) nebular line emission (panel a), but the most striking 
feature is the high He n /14686 / H/3 = 0.03 ratio (left, panel b 
and top ima ge). The FORS2 sp ectrum of C2_35, classified as an 
early OI bv ICastro etafl J2008h . is similar to Cl_31. The broad 
He H/14686 wind feature is indicative of the presence of a hot 
WR star. On top of the broad wind feature, there is a narrow 
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Figure 8. Upper panel: the predicted T(Om) from the synthesized neb- 
ular spectrum as a fu nction of T^g and L„ of a central ionizing O-star 
<Lanz & Hubenvl2003l) in a CLOUDY model with an inner radius of 0. 1 pc. 
The solid and dashed horizontal lines are the observed T(0 m) and error 
estimate for the direct environment of Cl_31. Lower panel: the predicted 
nebular line ratio He n A 4686/H/J as a function of r c fj and L, of a central 
ionising star. The black dots corresp ond to realistic dwarf and super giant 
O-type stars in the T^g and L* grid (Martins, Schaerer & Hillier 2005). 



nebular emission line. Hen emission lines from nebulae are only 
rar ely seen and often associated with strong X-ray sources (see 

e-g 
sec 



Pakull & Angebaulj 1 19861: iKaaret. Ward & Zezaj |2004) . but 



Shirazi & Brinchmann (2012). No obvious X-ray source is 



detected in archival XMM-Newton and Chandra observations at 
the location of C2_35 (priv. comm. R. Wijnands). The exposure 
times of these images, however, would not be sufficient to detect 
the X-ray emission of, for example, a stellar mass black hole at 
this distance. 



5 MODEL H II REGION 

If the stellar spectrum is a composition of different sources, various 
solutions are possible. However, the properties T(Om), n c , and 
[O/H] of the surrounding nebula are constrained (see Section [4]l, 
and the nebular emission profiles along the slit give some idea of 
the size of the ionised region. The ionising source, of which the 
constituents are constrained by the stellar spectrum, should be able 
to produce a region with properties we derive from the nebular 
spectrum. In this section we will use the s pectral synthesis code 
CLOUDY (version 08.00. iFerland et al.ll 19981) to investigate which 
properties of the nebula and of the central ionising source have a 
strong influence on the observables that we measured. This will 
allow us to constrain T e g and L, of the ionising source. 

CLOUDY is designed to simulate gaseous interstellar media. 
From a given set of conditions such as luminosity and spectral 
shape of the ionizing source, the program computes the thermal, 
ionisation and chemical structure of a region as well as the 
emitted spectrum. In order to compare our model results to what 



has been observed, we will mainly use the simulated emitted 
spectrum. Since we only know the total abundance of oxygen, 
we use the abundance pa ttern corresponding to the Orion Nebula 
provided by CLOUDY ^Baldwin et all Il99ll; iRubin et alj Il99ll; 



lOsterbrock Tran & Veilleuxll 19921; ISavage & Sembachlll996l) . and 
scale the metallicity such that the oxygen abundance matches our 
observed value. The total hydrogen density «h is set to 20 cirr 3 , 
consistent with the low electron density limit derived from [On] 
and [S n], assuming that all hydrogen is ionised. We do not include 
dust grain^ We use a spherical geometry: the inner radius of the 
cloud is 0. 1 pc and the outer radius is set where the temperature 
drops below 4 000 K. 

We examined the ratio R (Eq. [2]l of the modeled output 
spectra, and computed r(Ora), as a function of 7" eff and L, 
of a synthetic O-star spectrum dLanz & Hubenvl l2003h as cen- 
tral ionising source; see the upper panel of Fig. [8] A grid of 
T cff = 30 000 - 55 000 K and log(L/L Q ) = 3 - 6 has been examined; 
the black dots indicate realistic stars in this grid according to the 
calibration by iMartins, Schaerer & HillieJ d2005l) . We see that a 
higher r eff of the central star leads to a higher T(0 in), while the 
effect of L, on the 'measured' T(0 m) is smaller. L, does however 
influence the size of the cloud: a higher L, leads to a larger cloud. 
A cloud with metallicity Z = 0.3 Z needs a r cff ~ 44 000 - 55 000 
K star, depending on the luminosity, to produce the measured 
r(Oin) of 11 500 ± 600 K (horizontal line). Furthermore, we find 
that the metallicity of the cloud has an even stronger influence on 
T(Om); lower metallicities result in hotter clouds because they 
are less efficiently cooled. Changing «h only influences the spatial 
scale of the cloud: a density ten times as low results in a region 
four times as large. 

We also analysed the nebular line ratio He n A 4686/H/? as a 
function of r efl and L, of a central ionising source; see the lower 
panel of VigM Below T cff =* 40 000 K, no significant Hen A 4686 
line is predicted for any of the luminosities in our grid. For r cff > 
40 000 K, the line can be produced, and is stronger with respect to 
HyS for more luminous sources. The ratio He n A 4686/H/? decreases 
strongly by increasing the inner radius of the model cloud. The 
nebula directly around Cl_31 does not show Heii/i4686, therefore 
the inner radius is likely to be 1 pc or more. T(0 m) is not affected 
significantly by changing the inner radius. 



6 DISCUSSION: A CONSISTENT PICTURE 

In the previous sections, we have put constraints on the stellar 
parameters of Cl_31 and its surrounding region. We suggest that 
Cl_31 is not a single object, but rather a stellar cluster. We sum- 
marise the main arguments below. 

(a) the line profiles can not be reproduced by one single stellar 
atmosphere model, especially not Heii/14686 (Section[37L4}; 

(b) the visual absolute magnitude M v = -8.7 of this source is 
very high for a single object (Section[3}; 

(c) a very hot central object (~50000K) is necessary to pro- 
duce a T(0 m) of ~ 1 1 500 K in the surrounding nebula (SectionlS}, 
but the "average" spectral type suggests that the majority of the 



3 We exclude dust grains to keep the model simple. The extinction we mea- 
sure both in the stellar spectrum and the hydrogen emission lines could as 
well be due to dust that is outside the surrounding ionised region. 
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Figure 9. The observed normalised spectrum (solid lines) together with the 
combined line profiles (dotted lines) from the models described in Table [4] 
representing the late O giant component MOD89 (dash-dot-dot lines) and 
the WN-like component MOD90 (dashed lines), in the ratio 4:1. 



luminosity in the visual is produced by r c(f < 35 000 K stars (Sec- 
tion l3TTV 

The spectrum of a composite object results in a superposition 
of all components, weighed by their relative brightness at the 
wavelength considered. All lines that we analysed are close in 
wavelength, so we adopt a general ratio in brightness which is 



„ 10" 

v 10 1/ 
E 
o 

i 10 1t 

CL 

i 10 s 



- — T'i 



Lyman continuum 



input spectnjm 1 0'MOD89 
- - input spectnjm 2'MOD90 
total input spectrum 1 0*MOD89+2*MOD90 



100 



1000 



X (nm) 



Figure 10. The spectral energy distribution at the inner radius of the model 
cloud (r m = 0.1 pc), split in the different components. It is clear that in 
the visual, the 30 000 K component (MOD89, dot-dashed lines) dominates 
the light, while the hydrogen ionising photons (left of the vertical dashed 
line, with hv > 13.6eV) are provided mainly by the hot (50000 K) WN 
component (MOD90, dashed lines). The dotted profile shows the total input 
spectrum. 



the ratio in V. Given a universal initial mass function (IMF), a 
cluster will consist of many cool low-mass stars, and only a few 
very luminous and hot ones. In this analysis we will focus on the 
most massive and luminous members, because these dominate the 
cluster spectrum, as well as the flux of ionising photons. To this 
end we present a simple combination of FASTWIND models that 
(a) reproduces all observed line profiles, (b) has an absolute mag- 
nitude My = -8.7, and (c) when put into a CLOUDY model of an 
Hn region, produces an electron temperature r(Om) ~ 11 500 K, 
for an adopted metallicity Z = 0.3 Z Q and hydrogen density 
«h = 20 cirT 3 . 

Model MOD90 (see Table |4) has a high temperature 
(r eff = 50000 K), a high mass-loss rate (3 x 10- 5 M o yr-») 
and an enhanced helium abund ance Ny^INw = 0.8. It mimics a 
Wolf-Rayet WN star (see e.g. ICrowthe3 r2008). These properties 
result in a strong and broad He n /14686 emission feature (Fig. [9] 
upper right). To reproduce the observed shape of Heii/14686, 
we combine this profile with a model with a weak Heii/14686 
absorption profile: model MOD89, with T eff = 30000 K. This 
model resembles a late-O/early-B giant or bright giant. We create a 
combined profile of 20% MOD90 and 80% MOD89. With this flux 
ratio, the He i and He n absorption lines resemble a 30 000 K star, 
because in this respect the MOD89 model is dominant. The Ha 
wings are well reproduced (Fig. [9] upper right). The other B aimer 
lines in the WN component MOD90 are also affected by the strong 
wind, which results in a shallower line or even a P-Cygni profile in 
the case of H/3 and Hy. However, the strong absorption profile of 
the late O giant component MOD89 dominates, and the combined 
profiles match the observed ones. 

If the cluster consists of ten stars like MOD89 and two 
like MOD90, the ensemble would have a visual magni- 
tude of My = -8.63, in agreement with the observed value 
My = -8.7 ± 0.4. The flux ratio in the visual would be ~ 4 : 1, due 
to the different bolometric corrections (see Table [6J. The visual 
flux is dominated by the late O-type component, while further to 
the UV, the hot Wolf-Rayet component would dominate. The latter 
is required to reproduce the observed electron temperature in de 
cloud. 

We have used the combined spectrum of 10 times MOD89 
and 2 times MOD90 as ionising source in a CLOUDY model. 
The SED at the inner radius of the cloud is shown in Fig. [10] 
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Table 6. Properties of the final composition of the cluster. Our model clus- 
ter has two different components: a late O giant component (MOD89) and a 
Wolf Rayet WN-type component (MOD90), see also Tableg] Model ID, ef- 
fective temperature, bolometric luminosity, bolometric correction, absolute 
visual magnitude, number of stars of this type in the cluster, total absolute 
magnitude of this component, and fraction of the total visual flux provided 
by this component. 



ID 


Teff 


U 


BC 


My 


# 


My 


Fraction 




(K) 


(Lo) 








total 


of F vis 


MOD89 


30000 


10 5.24 


-2.42 


-5.93 


10 


-8.43 


0.8 


MOD90 


50000 


10 5.83 


-3.90 


-5.95 


2 


-6.70 


0.2 



where we see that the ionising flux is indeed dominated by the 
WN component MOD90. We use the same configuration as we 
did in Section [5] with Z = 0.3 Z Q . We choose n H = 20 crrT 3 , 
such that CLOUDY produces an ionised region with a radius of 
~20pc, see Section |4~T1 From the synthesized nebular spectrum 
we infer r(Om) ~ 10 800 K, slightly lower than the measured 
T(0 in) ~ 1 1 500 ± 600 K, but in reasonable agreement. According 
to the model in Section|3]and its results in Figure(8] a 50 000 K star 
would be able to heat the cloud to ~ 1 1 500 K. However, adding 
more late O stars decreases r(Om). The predicted nebular line 
ratio He n A 4686/H/J for this cluster composition is 0.06, which 
challenges the non-detection of nebular line Hen A 4686 around 
CT_31. This disagreement can be reconciled by increasing the 
inner radius of the model cloud to 1 pc or more. 

In principle, information as to the stellar content may also be 
derived from considering the mass-loss rates of the contributing 
stars. For our late OII/III source (MOD89) the adopted mass-loss 
rate of 2 x 10~ 6 M yr~' is rather large compared to theoretical ex- 
pectations for such a star at a metallicit y of 0.3 Z , being an order 
of ma gnitude higher than predicted by IVink. de Koter & Lamerj 
after co rrecting the empirica l mass-loss rate for wind in- 
homogeneities jMokiem et alj [2007h . Using the Vink et al. pre- 
scription for the WN star (MOD90) yields a much smaller dis- 
crepancy of a factor ~2. The discrepancy in the O star mass- 
loss rate can be partly reconciled by taking a smaller number of 
brighter stars, for instance supergiants, in the following way: A 
brighter star has a larger radius. In order to preserve the Ha pro- 
file shape, the quantity Q oc A<f/( j?^ 2 v 00 ) needs to remain invariant 
dde KoterTTl eap & Hubeny 1998). This implies that for fixed tem- 
perature M oc £,» , However, the expected mass-loss rate scales as 
M oc Lj- 2 , therefore a smaller number of brighter O stars may still 
match the strong Ho- emission line wings and better reconcile ob- 
served and theoretical mass-loss rates. We did not pursue this strat- 
egy in view of the uncertainties that are involved, for instance those 
relating to the stellar mass (which enters the problem as mass-loss 
is expected to scale with mass as M oc M~ 13 ), Moreover, we re- 
mark that higher than expected mass-l oss rates have been r eported 
for O stars in low-metallicity galaxies dTramper et alj|20llh . 



7 SUMMARY AND CONCLUSIONS 

We have analysed the VLT/X-shooter spectrum of Cl_31, one 
of the most luminous sources in NGC 55, and its surroundings. 
We conclude that NGC 55 Cl_31 is a cluster consisting of several 
massive stars, including at least one WN star, of which we observe 



the integrated spectrum. 

The H, He i and Hen lines in the stellar spectrum have 
been compared to synthesized spectra from a grid of FAST- 
WIND non-LTE stellar atmosphere models. All normalised 
lines except Hen A 4686 can be reproduced by a single-star 
model with T cff < 35 000 K, M ~ 2 x 10" 6 M Q yr 1 , and 
v rot sin(i) = 150 ± 50 fans" 1 . Heii/14686 has an equivalent width 
of -3.6 ± 0.4 A, but is ~3000km s _1 wide. No single star model is 
able to produce matching profiles for all lines simultaneously. 

Analysis of the nebular emission spectrum along the 
slit yields an electron density n c < 10 2 cm 4 , electron tem- 
perature T (O in) = 11 500 ± 600 K, and oxygen abundance 
[O/H] = 8.18 ± 0.03, which corresponds to a metallicity 
Z = 0.3 1 ± 0.04 Z Q . A grid of CLOUDY models suggests that a hot 
(~50 000 K) ionising source is necessary to reproduce the observed 
T (O in) in a H n region with comparable density and metallicity. 

We have also presented an illustrative cluster composition that 
reproduces all observed spectral features, the visual brightness of 
the target, and which is able to maintain an H n region with prop- 
erties similar to those derived from the nebular spectrum. In our 
model, the cluster contains several blue (super)giants and one or 
more WN stars. While the proposed composition might not be 
unique, the presence of at least one very hot, helium rich star with 
a high mass-loss is a robust conclusion. High angular resolution 
imaging reaching a resolution of 0.05" (corresponding to a phys- 
ical distance of about 0.5 pc) would provide an improvement of a 
factor 10 to 20 compared to our seeing-limited observations and 
would help to constrain the composition of the cluster. This makes 
NGC55C1_31 a prime target for ELT-class telescopes combin- 
ing high angular resolution and integral field or multi-object spec- 
troscopy. 
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